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Abstract. The possible role of decoherence due to space-time foam is discussed within the 
context of two models, one based on string/brane theory, and the other based on properties of 
black hole horizons in general relativity. It is argued that the density matrix satisfies a dissipative 
master equation, primarily from the study of renormalization group flows in non-critical string 
theory. This interpretation of the zero mode of the Liouville field as time leads necessarily to 
the CPT operator being ill defined. One striking consequence is that the quantum mechanical 
correlations of pair states of neutral mesons produced in meson factories are changed from the 
usual EPR state. The magnitude of this departure from EPR correlations is characterised by 
a parameter uj. The predicted value of u) is very small or zero. However it is shown explicitly 
that the the non- vanishing of u is only a feature of the model based on string/brane theory. 



1. Introduction 

The important feature of General Relativity (GR), the classical field theory describing gravity, 
is the fact that space-time is not simply a frame of coordinates, on which events take place, but 
is itself a dynamical entity. The quantisation of the theory, poses a problem, given that the 
coordinates of space-time themselves appear "fuzzy" . The fuzziness of space-time is associated 
with microscopic quantum fluctuations of the metric field, which may be singular. For instance, 
one may have Planck size (10~^^ m) black holes, emerging from the quantum gravity "vacuum", 
which may give space-time a "foamy", topologically non-trivial structure. The semi-classical 
study of black holes in GR has led to the prediction of Hawking radiation. In classical GR simple 
solutions for black holes exhibit event horizons. This was interpreted quantum mechanically in 
terms of lack of communication between Hilbert spaces within and outside the horizon. This 
is at the basis of the information paradox[T]. Simple arguments relating to tracing over states 
then lead to a description in terms of non-pure density matrices. These observations have led to 
a vigorous debate concerning the validity of unitary evolution for quantum gravity. Initially S. 
Hawking believed that there was information loss; he and others have recently claimed [2] that 
information is not lost in the presence of a black hole, but rather it is entangled in a holographic 
way with the portion of space-time outside the horizon. There seems currently more and more 
adherents to this point of view. However in my opinion the situation is not resolved. This is 
primarily because we do not have a theory of quantum gravity. Consequently there are always 
certain loose ends in the arguments that are proposed. As an example the original argument 
of Hawking has been criticized [3] on the grounds that horizons are mathematical constructs 
and are not as physically relevant as trapping surfaces. It is then argued that evaporating 
black holes can evade the information paradox owing to the presence of an inner and outer 
trapping surface. The recent arguments of Hawking and other (string) theorists are somewhat 



special in their details and require constructions such as the extremality of black holes, anti-de 
Sitter (supersymmetric) space-times and the Euclidean formulation for summing over space- 
time geometries. Such considerations have led to a point of view labelled as holography and 
black hole complementarity [3], [5] .The special features in these arguments make it hard to make 
general statements about unitarity or the lack of it in quantum evolution. 

We will take a somewhat different point of view which will put us firmly in the camp of those 
who believe that there is non-unitary evolution due to quantum gravity. Much of the discussion 
using string theory is based on critical strings where the demands of world sheet conformal 
and target space Lorentz invariance constrain the allowed values of space-time dimension (see 
e.g. [B]). Consequently changes in space-time backgrounds cannot be accommodated. Surely a 
proper theory of quantum gravity needs to be able to handle this. This is a fundamental problem. 
An approach which takes some steps to address these issues is based on the theory of non- 
critical strings [7] also referred to as two dimensional gravity coupled to matter. The restoration 
of conformal invariance for a number of spatial dimensions different from the critical value is 
addressed through the introduction of the Liouville field. Furthermore by studying the recoil of 
branes and the induced backreaction on space-time, it possible to consider some situations where 
the space-time metric and hence geometry is changed. When p space dimensional Dirichlet (Dp) 
branes are considered, it is no longer necessary to concentrate on the properties of microscopic 
black holes and the inherent issues of trapped surfaces. In particular Dp branes with p = 
exist in some string theories such as bosonic, type IIA and type I. Furthermore even when 
elementary D particles cannot exist consistently there can be effective D-particles formed by 
the compactification of higher dimensional D branes. Moreover D particles are non-perturbative 
constructions since their masses are inversely proportional to the the string coupling Qs ■ The 
introduction of the Liouville field opens up also an interesting possibility. It was suggested by 
Shore [8j that a renormalization group scale which was local on the world sheet was a useful 
way of analyzing string theory sigma models (in the sense of the Zamolodchikov C-function) and 
subsequently [7] this scale was identified with the Liouville field. Furthermore the zero mode 
of the Liouville field can be identified with time on the basis of transformation properties of 
vertex operators associated with D particles [9]. The Zamolodchikov C theorem (for general two 
dimensional field theories) requires that the C-function (a function in the space of all couplings, 
i.e. theory space) is a non-negative and non-increasing function on a renormalization group 
trajectory. The renormalization group trajectory can be identified with temporal evolution as 
just intimated and this leads, because of the increase of entropy along the trajectories, to a 
non-unitary evolution for the density matrix p. It will suffice for us to note that the evolution 
has the form 



Here Gij is the so called Zamolodchikov metric, /3* is the renormalization group /3 function 
for the couplings gi associated with vertex functions pj spanning theory space. Although the 
identification of the Liouville field with the local renormalization group scale is not rigorous, 
it is one important feature which influences our joining of the "non-unitary" camp. The other 
important reason for us is D-particle recoil which will be addressed in a later section. 

The second model to be considered by us is concerned with the effect of horizons and the 
consequent absence of unitarity but the formulation is not supported by a formal theory like 
string theory. A different effective theory of space-time foam has been proposed by Garay |10j . 






The fuzziness of space-time at the Planck scale is described by a non-fluctuating background 
which is supplemented by non-local interactions. The latter reflects the fact that at Planck scales 
space-time points lose their meaning and so these fluctuations present themselves in the non- 
fluctuating coarse grained background as non-local interactions. These non-local interactions 
are then rephrased as a quantum thermal bath with a Planckian temperature. The quantum 
entanglement of the gravitational bath and the two meson (entangled) state is explicit in this 
model. Consequently issues of back reaction can be readily examined. Since the evolution 
resulting from the standard Lindblad formulation does not lead to the uj effect, this manifestation 
of CPTV is not the result of some arbitrary non-unitary evolution. Hence it is interesting to 
study the above two quite distinct models (one motivated by string theory and the other by the 
properties of black holes in general relativity) for clues concerning the appearance of (and an 
estimate for the order of magnitude) of to. 

Recently these fundamental issues have been brought into sharp focus through a 
phenomenological analysis involving correlated neutral meson pairs. For equations such as ([T|), 
it is well known that initial pure states evolve into mixed ones and so the S-matrix jS relating 
initial and final density matrices does not factorise, i.e. 

^ / SS^ (3) 

where S = e^^^ . This is another way of saying that we have non-unitary evolution. In these 
circumstances Wald [TT] has shown that CPT is violated, at least in its strong form, i.e. there 
is no unitary invertible operator B such that 
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This result is due to the entanglement of the gravitational fluctuations with the matter system. 

It was pointed out in [12j, that if the CPT operator is not well defined then this has 
implications for the symmetry structure of the initial entangled state of two neutral mesons 
in meson factories. Indeed, if CPT can be defined as a quantum mechanical operator, then the 
decay of a (generic) meson with quantum numbers J^^ = 1 [13], leads to a pair state of 
neutral mesons |i) having the form of the entangled state 
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This state has the Bose symmetry associated with particle-antiparticle indistinguishability 
CV = +, where C is the charge conjugation and V is the permutation operation. If, however, 
CPT is not a good symmetry (i.e. ill-defined due to space-time foam), then Mq and Mq may not 
be identified and the requirement of CV = + is relaxed [I2J. Consequently, in a perturbative 
framework, the state of the meson pair can be parametrised to have the following form: 
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and Lo = \uj\ e'^ is a complex CPT violating (CPTV) parameter. For definiteness in what 
follows we shall term this quantum-gravity effect in the initial state as the "w-effect" |12j . There 
is actually another dynamical "w-effect" which is generated during the time evolution of the 
meson pair but this will not be discussed here [Hj. 

The structure of the article will be the following: we commence our analysis by discussing 
the omega effect. We will then proceed to discuss a string based model for the omega effect. The 
difficulty of generating the omega effect using other plausible approaches incorporating open 
systems and non-unitary master equations is illustrated through a class of models which we dub 
the thermal bath model. 



2. D- particles 

String theory was found, in the so called first revolution, to contain quantum gravity while in 
the second revolution the proposal of D(irichlet) brane solitons [151 116j . together with dualities 
brought about relationships between the different string theories that had been proposed earlier. 
In particular zero dimensional D-branes [T7] occur (in bosonic and some super symmetric string 
theories) and are known as D-particles. The spectrum of open strings attached to a Dp 
brane (with p > ) contains a Maxwell field and the ends of the open string carry charge. 
The associated Maxwell field is confined to the world volume of the D brane. Hence a 
conventionally charged string cannot end on a D-particle. We will thus restrict our consideration 
to neutral particles that are "captured" by D-particles. The D-particle solitons are fundamental 
scalar particles in the dual theories. Interactions in string theory are, as yet, not treated 
as systematically as in ordinary quantum field theory where a second quantised formalism is 
defined. The latter leads in a systematic way to the standard formulations by Schwinger and 
Feynman of perturbation series. When we consider stringy matter interacting with other matter 
or D-particles, the world lines traced out by point particles are replaced by two-dimensional 
world sheets. World sheets are the parameter space of the first quantised operators ( fermionic 
or bosonic) representing strings. In this way the first quantised string is represented by actually 
a two dimensional (world-sheet) quantum field theory. An important symmetry of this first 
quantised string theory is conformal invariance and the requirement of the latter does determine 
the space-time dimension and/or structure. This symmetry leads to a scaling of the metric and 
permits the representation of interactions through the construction of measures on inequivalent 
Riemann surfaces [IB] . In and out states of stringy matter are represented by punctures at the 
boundaries. The D-particles as solitonic states ||I6| in string theory do fiuctuate themselves, 
and this is described by stringy excitations, corresponding to open strings with their ends 
attached to the D-particles. In the first quantised (world-sheet) language, such fiuctuations 
are also described by Riemann surfaces of higher topology with appropriate Dirichlet boundary 
conditions (c.f. fig.[l]). The plethora of Feynman diagrams in higher order quantum field theory 
is replaced by a small set of world sheet diagrams classified by moduli which need to be summed 
or integrated over [6|. In order to understand possible consequences for CPT due to space-time 
foam we will have to characterise the latter. The model we will consider is based on D-particles 
populating a bulk geometry between parallel D-brane worlds. The model is termed D-foam |19] 
(c.f. figure E]), since our world is modelled as a three-brane moving in the bulk geometry. As 
a result, D-particles cross the brane world, and thereby appear as foamy flashing on and off 
structures for an observer on the brane. 

Even at low energies E, such a foam may have observable consequences e.g. decoherence 



or change in the usual Lorentz invariant dispersion relations. The study of D-brane dynamics 
has been made possible by Polchinski's realisation that such solitonic string backgrounds can be 
described in a conformally invariant way in terms of world sheets with boundaries [16]. On these 
boundaries Dirichlet boundary conditions for the collective target-space coordinates of the soliton 
are imposed [20.J. Heuristically, when low energy matter given by a closed string propagating in 
a (d + l)-dimensional space-time collides with a very massive D-particle (0-brane) embedded in 
this space-time, the D-particle recoils as a result [21] ■ We shall consider the simple case of bosonic 
stringy matter coupling to D-particles and so we can only discuss matters of principle and ignore 
issues of stability. However we should note that an open string model needs to incorporate for 
completeness, higher dimensional D-branes such as the D3 brane. This is due to the vectorial 
charge carried by the string owing to the Kalb-Ramond field. Higher dimensional D-branes 
(unlike D-particles) can carry the charge from the endpoints of open strings that are attached 
to them. For a closed bosonic string model the inclusion of such D-branes is not imperative (see 
figure [2]). The higher dimensional branes are not pertinent to our analysis however. The current 
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where Mp is the Planck mass 
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Figure 1. Upper picture: A Fluctuating D-particle is described by open strings attached to 
it. As a result of conservation of fluxes [151 HH [IT] that accompany the D-branes, an isolated 
D-particle cannot occur, but it has to be connected to a D-brane world through flux strings. 
Lower picture: World-sheet diagrams with annulus topologies, describing the fluctuations of 
D-particles as a result of the open string states ending on them. Conformal invariance implies 
that pinched surfaces, with infinitely long thin strips, have to be taken into account. In bosonic 
string theory, such surfaces can be resummed \22\ . 
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Figure 2. Schematic representation of a D-foam. The figure indicates also the capture/recoil 
process of a string state by a D-particle defect for closed (left) and open (right) string states, 
in the presence of D-brane world. The presence of a D-brane is essential due to gauge flux 
conservation, since an isolated D-particle cannot exist. The intermediate composite state at 
t = 0, which has a life time within the stringy uncertainty time interval 6t, of the order of the 
string length, and is described by world-sheet logarithmic conformal field theory, is responsible 
for the distortion of the surrounding space time during the scattering, and subsequently leads to 
induced metrics depending on both coordinates and momenta of the string state. This results 
on modified dispersion relations for the open string propagation in such a situation [19], leading 
to non-trivial "optics" for this space time. 



state of phenomenolgical modelling of the interactions of D-particle foam with stringy matter 
will be briefly summarised now. Since there are no rigid bodies in general relativity the recoil 
fluctuations of the brane and their effective stochastic back-reaction on space-time cannot be 
neglected. As we will discuss, D-particle recoil in the "tree approximation" i.e. in lowest order 



in the string coupling g^^ corresponds to the punctured disc or Riemann sphere approximation in 
open or closed string theory, induces a non-trivial space-time metric. For closed strings colliding 
with a heavy (non-relativistic) D-particle the metric has the form [23) 

9ij = Sij, 900 = -l,50i = e + Uit) @eit), i = l,...,d. (6) 

Here the suffix denotes temporal (Liouville) components, ^ oc 9s^^, where Aki = {ki — k2)^ 
and ki (^2) is the momentum of the propagating closed-string state before (after) the recoil, 
fji are the spatial collective coordinates of the D- particle and is identified with the target 
Minkowski time t for t S> after the collision; {t) is the regularised step function represented 
by 

1 r dq 



ee{t) = — —^e^'\ (7) 
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and Ui is small,. For our purposes the Liouville and Minkowski times can be identified. Now for 
large t, to leading order, 

90i " = Y 

where Apj is the momentum transfer during a collision and Mg is the string mass scale, (7<j < 1 
is the string coupling, assumed weak, and the combination Mg/gs is the D-particle mass, 
playing the role of the Quantum Gravity scale in this problem, i.e. the Planck mass; this 
formalism has been used to establish a phenomenological model where the couplings Ui are 
taken to be stochastic and modeled by a gaussian process[24j. The stochasticity is due to two 
contributions. The D-particle has an initial random recoil velocity due to capture of stringy 
matter and subsequent release. Furthermore sumperimposed on this randomness are quantum 
fluctuations due to vacuum string excitations. The process of capture and emission does not have 
to conserve flavour. Consequently we need to generalise the stochastic structure to allow for 
this. The fluctuations of each component of the metric tensor g'"^ will then not be just given by 
the simple recoil distortion ([61), but instead will be taken to have a 2 x 2 ("flavour") structure 

= (-1 + ^-4)1 

5°^ = 5^° = '^oi + riCTi + r2(T2 + r^a^ 

9'' = (l+rs)! (9) 

where 1 , is the identity operator and cjj are the Pauli matrices and the are random variables 
which obey 

(r^) = 0, (r^r^) = A^J^^. (10) 

The above parametrisation has been taken for simplicity and also we will consider motion to 
be in the x- direction which is natural since the meson pairs move collinearly. The above 
parametrisation has been taken for simplicity and we will also consider motion to be in the x- 
direction which is natural since the meson pairs move collinearly. In any given realisation of 
the random variables for a gravitational background, the system evolution can be considered 
to be given by the Klein-Gordon equation for a two component spinless neutral meson field 

= ^ 0^ ^ (corresponding to the two fiavours) with mass matrix m = \ (mi + m2) 1 + 
^ (mi — m2) (T3. We thus have 

{g'^^D^Dp - m2)$ = (11) 

where is the covariant derivative. Since the Christoffel symbols vanish (as the ai are 
independent of space-time) the coincide with da- Hence within this fiavour changing 



background (llip becomes 
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We shall consider the two particle tensor product states made from the single particle states 
t'^*^^ and 1^*^^ i = 1,2. This includes the correlated states which include the state of the uj 
effect. In view of (|12p the evolution of this state is governed by a hamiltonian H 
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which is the natural generalisation of the standard Klein Gordon hamiltonian in a one particle 
situation where k \ ±k, a) = zizk \ ±k, a) with a =| or |. H is a single particle hamiltonian and in 
order to study two particle states associated with i = 1, 2 we can define Hi in the natural way (in 
terms of H) and then the total hamiltonian is = X]i=i ^i- The effect of space-time foam on 
the initial entangled state of two neutral mesons is conceptually difficult to isolate, given that the 
meson state is itself entangled with the bath. Nevertheless, in the context of our specific model, 
which is written as a stochastic hamiltonian, one can estimate the order of the associated w-effect 
of [H] by applying non-degenerate perturbation theory to the states T^*^^, i = 1, 2 (where 

the label ±k is redundant since i already determines it). Although it would be more rigorous to 
consider the corresponding density matrices, traced over the unobserved gravitational degrees 
of freedom, in order to obtain estimates it will suffice formally to work with single-meson state 
vectors . 

Owing to the form of the hamiltonian (jl3p the gravitationally perturbed states will still be 
momentum eigenstates. The dominant features of a possible a;-effect can be seen from a term 
Hj in the interaction hamiltonian 
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(14) 



which is the leading order contribution in the small parameters rj (c.f. (I9[). ()13p ) in H (i.e yA 
are small). In first order non-degenerate perturbation theory the gravitational dressing of 
leads to a state: 
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and correspondingly for the dressed state is obtained from (fTB]) by \[) <-> |t) and a — > /? 

where 
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Here the quantities Ei = (m? + /c^)^/^ denote the energy eigenvalues, and i = 1 is associated 
with the up state and i = 2 with the down state. With this in mind the totally antisymmetric 
"gravitationally-dressed" state can be expressed in terms of the unperturbed single-particle 
states as: 
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It should be noted that for rj oc Sn the generated w-hke effect corresponds to the case 
q;(*) = — /3'^*\ while the w-effect of [25j ([5]) corresponds to oc 5i2 (and a^*) = /J^*^). The 
velocity recoil u is given by 

ti ~ r 

and so the variance A can be determined once the variance of u is known. In the perturbative 
approach we do not need to specify the higher order moments provided they are finite and small. 
Rigour would demand a density matrix approach but in order to estimate the order of the effect 
it suffices to note that 
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which is an estimate for |a;p. It is necessary to examine the moments for u. 



3. A string theory estimate for the variance of the D-particle 

In order to estimate the variance of the recoil velcoity we need to understand the dynamics of 
D-particle foam. It is imperative to consider D-particle quantisation and in particular D-particle 
recoil as a result of scattering off stringy matter. The issue of recoil is not fully understood and 
it is not our purpose here to delve into the subtleties of recoil and operators describing recoil. 
We will rather proceed on the basis of a proposal which has had some success in the past [21j. 
The energy of a D-particle is independent of its position. Consequently in bosonic string theory 
there are 25 zero modes, the Dirichlet directions. Zero modes lead to infrared divergences in 
loops in a field theory setting where collective co-ordinates are used to isolate these infrared 
divergences. This is essentially due to the naivety of the perturbation series that is used and can 
be addressed using a coherent state formalism. The situation is similar but, in some ways, worse 
for string theory since the second quantised formalism for strings is more rudimentary. In string 
theory the use of first quantisation requires a sum over Riemann sheets with different moduli 
parameters. The formal transition from one surface to another of lower genus (within string 
perturbation theory ) has singularities associated with infrared divergences in the integration 
over moduli parameters. In the case of a disc D an incipient annulus S can be found by making 
two punctures and attaching a long thin strip (somewhat like the strap in a handbag and dubbed 
wormholes). In the limit of vanishing width this wormhole, is represented by a region in moduli 
space, integration over which leads to an infrared divergence. A similar argument would apply 
to a Riemann sphere where two punctures would be connected by a long thin tube. For the disc 
the divergence will be due to propagation of zero mode open string state while for the Riemann 
sphere it would be an analogous closed string state along the wormhole. Let us examine this 
divergence explicitly. Consider the correlation function (V1V2V3 . . . ln)s for vertex operators 
[18j „ where S is the Riemann surface of an annulus; it can be deformed into a disc 

with a wormhole attached [20j.|22j. The set of vertex operators include necessarily any higher 
dimensional D-branes necessary to conserve string charge. However such vertex insertions clearly 



do not affect tlie infrared divergence caused by wormwholes. Tlie correlation function can be 
expressed as 



(^1^2 ^3 • • • K)s = E / / J ^/''-l {(l)a (si) (^2) 1^1 ^2 ^3 • • • K>^ (20) 

where si and S2 are the positions of the punctures on 9-D.The {(j)a} are a complete set of 
eigenstates of the Virasoro operator Lq with conformal weights /ia|26] and g is a Teichmuller 
parameter associated with the added thin strip. Clearly there is a potential divergence associated 
with its disappearance g — > and this corresponds to a long thin strip attached to the disc. For a 
static D-particle the string co-ordinates d = , , , . . . ,^^^} have Dirichlet boundary 
conditions while X^ has Neumann boundary conditions (in the static gauge) 
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where (r, cr) is a co-ordinisation of the worldsheet.The associated translational zero mode is 
given by 

f{X,co) = ^dnXhe'^''\ (22) 

where 9„ denotes a derivative in the X* Dirichlet direction, and is an element of the set {(pa}- 
The conformal weight hi = 1 + oiuP' . The relevant part of the integral in ()20p is 

1 /-co , 9 1 

, dq\ = / ^ (23) 

J-oo h g(-logg) ' 

and is divergent because of the behaviour of the integrand near q = ^. This can be regularized by 
putting a lower cut-off > 5 ^ in the integral. The correlation function on S can be computed 
to be 120! 



iyxV^Vz • • • ^n)E = log ^ {dn^D (si) .dn^ D (^2) ^1^21^3 • ■ ■ (24) 

where T is a cut-off for large (target) time. Division by it removes divergencies due to the 
integration over the world-sheet zero modes of the target time. These should not be confused 
with divergencies associated with pinched world-sheet surfaces, proportional to logd, that we 
are interested in here. These latter divergences cause conventional conformal invariance to fail. 

However, in one approach, it was argued sometime ago, that these divergences can be 
canceled if D-particles are allowed to recoil [2H [22], as a result of momentum conservation 
during their scattering with string states. In fact, in our D-particle foam model ^9j, this is 
not a simple scattering process, as it involves capture and re-emission of the string state by the 
D-particle defect. In simple terms, this process involves splitting of strings by the defect. In a 
world-sheet (first quantization) framework, such processes are described by appropriate vertex 
operators, whose operator product expansion close on a (local) logarithmic algebra |2HI22|. The 
translational zero modes are associated with infinitesimal translations in the Dirichlet directions. 
Given that D-particles do not have any internal or rotational degrees of freedom, these modes 
should give us valuable information concerning recoil. Moreover (for a; = 0) there is a degeneracy 
in the conformal weights between 4>'^, d^jcj)^ and the identity operator and also the conformal 
blocks in the corresponding algebra have logarithmic terms. 

To understand the formal structure of the world-sheet deformation operators pertinent to the 
recoil/capture process, we first notice that the world-sheet boundary operator Vd describing the 
excitations of a moving heavy DO-brane is given in the tree approximation by: 



^0= J {vidnX' + UiX'>dnX') = JVi (XO) dnX' 
dD dD 

where Ui and yj are the velocity and position of the D-particle respectively and Yi {X^) = 
Hi + UiX^ . To describe the capture/recoil we need an operator which has non-zero matrix 
elements between different states of the D-particle and is turned on "abruptly" in target time. 
One way of doing this is to put [21] a Q {X^), the Heavyside function, in front of Vd which 
models an impulse whereby the D-particle starts moving at X^ = 0. Using Gauss's theorem this 
impulsive Vd, denoted by Vd^^^, can be represented as 

25 „ 25 „ 

y^imp = Y^J^ d'^zda e d^x') = ^ j^^ druiX^e dnX\ (25) 

Since X^ is an operator it will be necessary to define {X^) as an operator using the contour 
integral 

G. (XO) = -f r ^e--° (26) 



with e — > 0+. 

Hence we can consider 

D,{X^) ^ D{X';e) = X^Q, (x^) = - T ^i^e-^° . (27) 

^ ^ J-oo {u> — ie) 

The introduction of the feature of impulse in the operator breaks conventional conformal 
symmetry, but a modified logarithmic conformal algebra holds. A generic logarithmic algebra 
in terms of operators C and T> and the stress tensor T (z) (in complex tensor notation ) satisfies 
the operator product expansion 

T (z) C («;, w) ~ -^C {w, w) + + ■■■ 

Tiz) D {w,w) ~ -^D {w,w) + -^C {w,w) + + ■■■ 

^ ' \ ' / (2:— «;) ^ ' (z—w) ^ ' yz—w) 

and 
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\z\ 



{D{z,z)D{Q,Q)) ~ ^{\og\z\+d), (30) 
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where d is a constant. Since the conformal dimension of e**^"^" is ^ we find that 

T {w) De (z) ~ - . . D, (z) + ^ . eOe M +■■■ (31) 

2[w — z) [w — z) ^ ^ 

and so a logarithmic conformal algebra structure arises if we define 

Ce{X^)^C{X^;e)=eQ,[x^) , (32) 



suppressing, for simplicity, the non-holomorphic piece. The above logarithmic conformal field 
theory structure is found with this identification. Similarly we find 

T{w)Ce{z)^- + 
2{w — z) 

2 

Consequently A for Ce {z) and (z) is — ^. A calculation (in a euclidean metric) for a disc of 
size L with a short-distance worldsheet cut-off a reveals that as e — > 



{Ce{z,z)Ce{m) 
{Ce{z,z) D, (0,0)) 

{De iz,z) De (0,0)) 



vr 



vr 



2 V e^a 
vr 



-W^ l-4e^log 



vr 



1 



2 V e^a W 



4 log 



where a = log 



. We consider e — > 0-|- such that 



(33) 
(34) 

(35) 
(36) 



where ij is the time signature and the right-hand side is kept fixed as the cutoff runs; it is then 
straightforward to see that (|33]l . (f3i|l . and (f35]) are consistent with (f28l) . (f29l) . and (f30l) . It is 
only under the condition (f36jl that the recoil operators and obey a closed logarithmic 
conformal algebra |21j : 



(C,(z)C,(0)) ~ (37) 
{Csiz)D,{0)) - 1 (38) 



(Z)a^)^£(0)) ~ -477 log 



(39) 



The reader should notice that the full recoil operators, involving holomorphic pieces 

with the conformal-dimension-one entering (j25p ). obey the full logarithmic algebra (j28p . (j29p . 

2 

(|30p with conformal dimensions A = 1 — ^ . From now on we shall adopt the euclidean signature 
T] = 1. 

We next remark that, at tree level in the string perturbation sense, the stringy sigma model 
(inclusive of the D-particle boundary term and other vertex operators) is a two dimensional 
renormalizable quantum field theory; hence for generic couplings it is possible to see how 
the couplings run in the renormalization group sense with changes in the short distance cut- 
off through the beta functions In the world-sheet renormalization group [27], based on 
expansions in powers of the couplings, /3* has the form ( with no summation over the repeated 
indices) 

= Vig' + ... (40) 

where yi is the anomalous dimension, which is related to the conformal dimension Aj by 
Ui = Ai — 6, with 6 the engineering dimension (for the holomorphic parts of vertex operators for 
the open string 6 = 1). The ... in (|40p denote higher orders in g^. Consequently, in our case, we 
note that the (renormalised) D-particle recoil velocities constitute such cr-model couplings, 
and to lowest order in the renormalised coupling Ui the corresponding /3 function satisfies 



where A is a (covariant) world-sheet renormahzation-group scale. In our notation, we identify 

r 2 

the logarithm of this scale with a = log — , satisfying (j36p . 

An important comment is now in order concerning the interpretation of the flow of this world- 
sheet renormalization group scale as a target-time flow. The target time t is identified through 
t = 2 log A. For completeness we recapitulate the arguments of [21j leading to such a conclusion. 
Let one make a scale transformation on the size of the world-sheet 

L^L' = e*/^L (42) 

which is a finite-size scaling (the only one which has physical sense for the open string world- 
sheet). Because of the relation between e and L (|36p this transformation will induce a change 
in e 

2 ^ 12 ^ (43) 

(note that if e is infinitesimally small, so is e' for any finite t). From the scale dependence of the 
correlation functions (j39p that and Dg, transform as C^/ = and — > -D^' = D^+tCs- 

Hence the coupling constants in front of and in the recoil operator ([3]), i.e. the velocities 
Ui and spatial collective coordinates yi of the brane, must transform like: 

Ui^Ui , yi + Uit (44) 

This is nothing other than the Galilean transformation for the heavy D-particles and thus it 
demonstrates that the finite size scaling parameter t, entering (|42p . plays the role of target time, 
on account of ()36p . Notice that (I44p is derived upon using ()39p . that is in the limit where e — > 0. 
This will become important later on, where we shall discuss (stochastic) relaxation phenomena 
in our recoiling D-particle. 

Thus, in the presence of recoil a world-sheet scale transformation leads to an evolution of the 
D-brane in target space, and from now on we identify the world-sheet renormalization group 
scale with the target time t. In this sense, equation (f^T]) is an evolution equation in target time. 

However, this equation does not capture quantum-fluctuation aspects of about its classical 
trajectory with time Ui{t). Going to higher orders in perturbation theory of the quantum field 
theory at fixed genus does not qualitatively alter the situation in the sense that the equation 
remains deterministic. However the effect of string perturbation theory where higher genus 
surfaces are considered and re-summed in some appropriate limits can be shown to give rise to 
quantum fluctuations in [28] which leads to a Langevin equation [29j rather than (|4ip. The 
equation has the form 

^ = _|^* +^tV2^(t) (45) 
dt At .pM ' ^ ^ 



where t = e^^ and (t) represents white noise. This equation is valid for large t. From the above 
analysis it is known that [22] (c.f. Appendix A)) that to O (g'^) the correlation for (t) is 
independent, and for time scales of interest, is correlated like white noise ; hence the correlation 
of ^ (t) has the form: 

{at)Cit'))=6it-t'). (46) 

Since the vectorial nature of is not crucial for our analysis, we will suppress it and consider 
the single variable u. Moreover the mean of u vanishes. 



4. Langevin equations and fluctuating string coupling 

It is possible to solve |28j the Langevin equation (I45p (with (l46p ).The quantity of interest which 
is an important input for our estimate of the omega effect is the variance of u the velocity of the 
D-particle. If the D-particle is typically interacting with matter on a time scale of to, then the 
effect of a large number of such collisions can be calculated by performing an ensemble average 
over a distribution of uq the initial velocity at the time of capture; a gaussian distribution 
with zero mean and variance a is assumed for uq and is appropriate in the absence of a priori 
knowledge. Using generic properties of strings [30] consistent with the space-time uncertainties 
|31j . the capture and re-emission time to, involves the growth of a stretched string between the 
D-particle and the D-brane world (e.g. a D3 brane corresponding to our 3-dimensional world). 
The minimisation of the energy associated with intermediate string tension and oscillations leads 
to : 

to ~ a'p^ 

where is the energy of the intermediate string state. For neutral kaons it is appropriate to 
consider non-relativistic particles and so ~ m| -|- m^. This result is not modified even if 
there are more than 2 D-particles in close proximity. Using this we can infer [28j from the 
solution of the Fokker-Planck equation associated with Eqn. (j45| ) that 



var (u) ~ 



2(l + V2) [g^.a'pl + 15^7^] 



15 



for the variance of u due to capture of stringy matter by D particles. The terms of order ci^ have 
not been written explicitly for several reasons. In most models of quantum gravity, we expect 
that < g^a'po since it is a natural assumption to make for a dispersion due to (quantum) 
fluctuations of the recoil velocity of heavy D-particles of average mass Ms/gg.We now return to 
the estimate for the uj effect [191 It is straightforward to show that 



and 



Hence 



= k 



{4 



El 
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, ■ (2) 
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2r 
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on noting that 



To r. 




{1).{2) 



= 

i = 1,2. 



This embodies the independence of the random variables which differ in their indices. 
Furthermore at this order in perturbation theory no information is required about higher order 
moments (provided they remain bounded). From Eqn. (jl9p and Eqn.Q we can conclude that 



o(p|) 



fgk"^ m\ + ml 



2u2 



{mi - m2)^ Af| {mi - 7712)^ 



where ("^ = '^iHh^Uhl^ the estimate given in [T^ was undetermined. It is a success of this 
analysis that it is possible to obtain a simple formula for it. On remembering the restrictions 



oc 6ii or Tj (X 6i2 discussed earlier, we find for kaons tliat C = O (1) and ~ 2 x 10"^ . Strictly 
speaking, since the constituent quarks making up the kaons are charged, it is the gluons of the 
strong force which interact with the D particles. However we shall not pursue this issue here but 
assume, simplistically, that the kaon is an elementary neutral excitation. The present bound for 
from DA<I>NE is < 10"'^ with 95% confidence level. Hence the the strength of the effect 
is not far below the sensitivity of current experiments and may be detectable when DA$NE is 
upgraded. A similar analysis may be done for B mesons produced in Super B factories. 



5. The Thermal Master Equation 

Master equations with a thermal bath have been argued to be relevant to decoherence with space 
time foam jlOj . His argument was based on problems of measurement [36]. Also, by considering 
the infinite redhshifting near the horizon for an observer far away from the horizon of a black 
hole, Padmanabhan [37] has argued that a foam consisting of virtual black holes would magnify 
Planck scale physics for observers asymptotically far from the horizon and thus an effective 
non-local field theory description should emerge. Garay has argued that mildest non-locality, 
bilocality, can be rewritten with the help of a local auxilary field and this plausibly leads to a 
description of space-time foam as a heat bath. A thermal field represents a bath about which 
there is minimal information since only the mean energy of the bath is known, a situation valid 
also for space-time foam. In applications of quantum information it has been shown that a 
system of two qubits (or two-level systems) initially in a separable state (and interacting with a 
thermal bath) can actually be entangled by such a single mode bath [38]. As the system evolves 
the degree of entanglement is sensitive to the initial state. The close analogy between two-level 
systems and neutral meson systems, together with the modelling by a phenomenological thermal 
bath of space time foam, makes the study of thermal master equations a rather intriguing one 
for the generation of to. The hamiltonian 7i representing the interaction of two such two level 
atoms with a single mode thermal field is 



aaV + a'a_ 



i=l 



(47) 



where a is the bosonic annihilation operator for the mode of the thermal field and the a's are 
the standard Pauli matrices for the 2 level systems. The operators a and satisfy 



a, a 



1, a\a^ 



[a,a\ 



0. 



(48) 



The superscripts label the particle. This hamiltonian is known as the Jaynes-Cummings 
hamiltonian [39j and explicitly incorporates the back reaction or entanglement between system 
and bath. This is in contrast with the Lindlblad model and the Liouville stochastic metric model. 
In common with the Lindblad model it is non-geometric. In the former the entanglement with 
the bath has been in principle integrated over while in the latter it is represented by a stochastic 
effect. An important feature of TC in[37]is the block structure of subspaces that are left invariant 
by TC. It is straightforward to show that the family of invariant irreducible spaces £n may be 
defined by | Cj-"^^ ,i = 1, . . . ,4| where ( in obvious notation, with n denoting the number of 
oscillator quanta) 



(49) 



and 



)s|l<')a<«.n + 2) 



The total space of states is a direct sum of the £n for different n. We will write 7i as + 
where 



and 



(50) 
(51) 



1=1 



n is a quantum number and gives the effect of the random environment. In our era the strength 
7 of the coupling with the bath is weak. We expect heavy gravitational degrees of freedom and 
so ^ V. It is possible to associate both thermal and highly non-classical density matrices 
with the bath state. Nonethelss because of the block structure it can be shown completely 
non-perturbatively that the uj contribution to the density matrix is absent |l0]. 

We will calculate the stationary states in using degenerate perturbation theory, where 
appropriate. We will be primarily interested in the dressing of the degenerate states 



and 



in) 



An) 



because it is these which contain the neutral meson entangled state . In 2nd order 



perturbation theory the dressed states are 
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+ 7' 
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with energy E. 
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(n + 1) 1/ + O (7'3) and 
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(52) 



(53) 



4"^ 



with energy E'^'' = {n + I) + ^ + O (7'^). It is 

effect. More precisely we would construct the state Tr ( (^2"^ 
density matrix (and has the form ps = J2nmPnn' \n) (n'| for suitab 



which can in principle give the iv- 
Pb^ where pb is the bath 
e choices of pnn' ) • To this 



and 



An) 



order of approximation 7/^2"^ ) cannot generate the w-effect since there is no admixture of 



An) 



However it is a priori possible that this may change when higher orders in 7. We 



can show that 



4"^ 



An) 



An) 



to all orders in 7 by directly considering the hamiltonian 



matrix for n within it is given by 
T-Cin) 



( + ni/ 7\/^ + 1 7\/?T- + 1 
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7Vri+ 1 (n^\)v 
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We immediately notice that 
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(55) 



and so, to all orders in 7, the environment does not dress the state of interest to give the w-effect; 
clearly this is independent of any choice oi ps- 



As we noted this block structure of the Hilbert space is a consequence of the structure of 
Til which is commonly called the hamiltonian in the 'rotating wave' approximation [IT] . (This 
nomenclature arises within the context of quantum optics where this model is used extensively.) 
It is natural to expect that, once the rotating wave approximation is abandoned, our conclusion 
would be modified. We shall examine whether this expectation materialises by adding io 7i a. 
'non-rotating wave' piece 7i2 



2 
i=l 



aa_ + a'aV 



(56) 



TC2 does not map £n into £n as can be seen from 



^2 



An) 
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,{n+2)\ 
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We have noted that 



annihilates 



An) 



An) 
An) 



A") 



is an eigenstate of TIq +Tli. However the perturbation Ti.2 



62 ' j and so this eigenstate does not get dressed by 7^2 • Hence the oj effect 
cannot be rescued by moving away from the rotating wave scenario. 

We should comment that a variant of the bath model based on non-bosonic operators a can 
be considered. This departs somewhat from the original philosophy of Garay but may have 
some legitimacy from analysis in M-theory. From an infinite momentum frame study it was 
suggested that D-particles may satisfy infinite statistics [42] . Of course these D-particles would 
be light contrary to the heavy D-particles that we have been considering and so we will consider 
a more general deformation of the usual statistics.. If fiavour changes can arise due to D-particle 
interactions then it may be pertinent to reinterpret the hamiltonian in this way. Generalised 
statistics is characterised conventionally by a c-number deformation [43] parameter q with g = 1 
giving bosons. Although we have as yet not made a general analysis of such deformed algebras, 
we will consider the following particular deformation of the usual harmonic oscillator algebra: 



'^a\aq 



(57) 



where the g-operators aq,a\^ and Nq act on a Hilbert space with a denumerable orthonormal 



basis , n = 0, 1, 2, . . .| as follows: 



4 =V^+ 1]^^^ 1^ + l)g ' «g = N^''^ \n - l)g > \n) q = n \n)^ (58) 



and 



(59) 
q-q ^ 

Clearly when q = 1 the boson case is recovered. We can consider the q generalisation of TC. The 
resulting q operator is 



Hq = nva\aq + hi^a'i^ + hi^af^ +^^1^2 («9^? + 4^- ) 



(60) 



1=1 



and there is still the natural q analogue Sn"^ of the subspace £n which is invariant under Tiq and 
spanned by the basis 

S"^) -|T(^\T(^))K,|er^) =|T«,i(^))|n + l)^,|e(")) ^ jj^, T^^)) |n + 1)^ , (61) 



and 



(62) 



A very similar calculation to the bosonic case yields 







Jy/[n + 1 











(63) 



Again the vector 



+ 1] 

4"^ ^ is an eigenvector with eigenvalue [n + 1] . Hence this model 

using generalised statistics based on the Jaynes-Cummings framework also does not lead to the 
w-effect. 

One cannot gainsay that other more complicated models of 'thermal' baths may display 
the w-effect but clearly a rather standard model rejects quite emphatically the possibility of 
such an effect. This by itself is very interesting. It shows that the w-effect is far from a generic 
possibility for space-time foams. Just as it is remarkable that various versions of the paradigmatic 
'thermal' bath cannot accommodate the w-effect, it also remarkable that the D-particle foam 
model manages to do so very simply. 



6. Conclusions 

In this work we have discussed two classes of space-time foam models, which may characterise 
realistic situations of the (still elusive) theory of quantum gravity. In one of the models, inspired 
by non-critical string theory, string mattter on a brane world interacts with D particles in the 
bulk. Recoil of the heavy D particles owing to interactions with the stringy matter produces 
a gravitational distortion which has a backreaction on the stringy matter. This distortion, 
consistent with a logarithmic conformal field theory algebra, depens on the recoil velocity of 
the D-particle. It is modelled by a stochastic metric and consequently can affect the different 
mass eigenstates of neutral mesons. There are two contributions to the stochasticity. One is due 
to the stochastic aspects of the recoil velocity found in low order loop amplitudes of open or 
closed strings ( performed explicitly in the bosonic string theory model) while the other is the 
leading behaviour of an infinite resummation of the subleading order infrared divergences in loop 
perturbation theory. The former is assumed to fluctuate randomly, with a dispersion which is 
viewed as a phenomenological parameter. The latter for long enough capture times of the D- 
particle with stringy matter gives a calculable dispersion which dominates the phenomenological 
dispersion for the usual expectations of the quantum gravity scale. This gives a prediction for the 
order of magnitude of the CPTV cj-like effect in the initial entangled state of two neutral mesons 
in a meson factory based on stationary (non degenerate) perturbation theory for the gravitational 
dressing of the correlated meson state. The Klein-Gordon hamiltonian in the induced stochastic 
gravitiational field was used for the calculation of the dressing. The order of magnitude of uj may 
not be far from the sensitivity of immediate future experimental facilities, such as a possible 
upgrade of the DA<I>NE detector or a B-meson factory. 

This causes a CPTV w-like effect in the initial entangled state of two neutral mesons in 
a meson factory, of the type conjectured in [12]. Using stationary perturbation theory it was 
possible to give an order of magnitude estimate of the effect: the latter is momentum dependent, 
and of an order which may not be far from the sensitivity of immediate future experimental 
facilities, such as a possible upgrade of the DA<I>NE detector or a B-meson factory. A similar 
effect, but with a sinusoidal time dependence, and hence experimentally disentanglable from the 
initial-state effect , is also generated in this model of foam by the evolution of the system. 



A second model of space-time foam, that of a thermal bath of gravitational degrees of 
freedom, is also considered in our work, which, however, does not lead to the generation of an 
w-effect. Compared to the initial treatment of thermal baths further significant generalisations 
have been considered here. In the initial treatment (using a model derived in the rotating 
wave approximation) the structure of the hamiltonian matrix was block diagonal, the blocks 
being invariant 4-dimensional subspaces.The real signifcance of the block diagonal structure 
was unclear for the absence of the omega effect in this model and has remained a matter for 
debate. We have analyzed the model ( now without the rotating wave approximation) and to 
our surprise the omega effect is still absent. Furthermore we have also considered q oscillators 
for the heat bath modes since from Matrix theory we know that energetic D-particles satisfy q 
statistics. For the model of q statistics that we have adopted there is again no w-effect in the q 
version of the thermal bath. 

It is interesting to continue the search for other (and more) realistic models of quantum 
gravity, either in the context of string theory or in alternative approaches, such as loop quantum 
gravity, exhibiting intrinsic CPT Violating effects in sensitive matter probes. Detailed analyses 
of global data in relation to CPTV, including very sensitive probes such as high energy neutrinos, 
is a good way forward. 
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